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CIRCULAR CYLINDRICAL SHEEi UNDER AXIAL COMPRESSION 
By Y. Yoshimura^ 


INTRODUCTION 


Ever since A. Robertson (ref. l), E. E. Lundqulst (ref. 2 .), and 
L. H. Donnell (ref. 3)^ shoved that the experimental buckling load of a 
circular cylindrical shell under axial con 5 )resslon is much smaller than 
the theoretical value much effort has been devoted tq this problem. 

Among the more recent theories, the nonlineEir finite deformation theory 
given by Th. von Kairman and Hsue-Shen Tsien (ref. 4) and T. Kono (ref. 5) 
explains convincingly the discrepancy between the experimental and theo- 
retical buckling load hitherto computed by the classical theory, and the 
large and abrupt displacement and the decrease of load in buckling. 

The stress-strain relationship obtained by Von K^man, etc., certainly 
indicates an important feature of the buckling of a circular cylindrical 
shell. However, the idea that the minimum value of the load given by the 
load-deflection curves may be realized is not considered appropriate.^ 

The state which may be actually realized after buckling must be deter- 
mined by minimizing the energy, not only with respect to the magnitude of 
deflection, but also to the aspect ratio and the circumferential niomber 
of buckling waves. The actual buckling load will be given by a comparison 
of energy levels before and after buckling and the energy barrier to be 
jumped over in buckling. 5 Based on such a concept, the general buckling 
and the local buckling of a cylindrical shell are considered to be quite 
different phenomena from the energy viewpoint, though they are equivalent 
with respect to the load. 


After this investigation was performed, the author became aware of 
Professor H. S. Tsien’ s paper "A Theory for the Buckling of Thin Shells", 
Journ. Aero. Scl . 9, 1942, 375-384, v/hich had remained unknown to us 
due to the war. 

2 

Associate Professor of University of Tokyo. 

3 

Note by reviewer: The energy level concept was also clarified by 

Von Karman and Tsien in the paper quoted above. 
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The reason why the cylindrical shell reveals such a particular 
behavior for buckling will be attributed to the geometrical features of 
its deformation, that is, to the fact that a developable surface quite 
different from the original cylindrical sijrface can exist. On the basis 
of these fundamental concepts concerning deformation and energy, the mech- 
anism of the buckling can be understood more completely. 


LIST OF SYMBOLS 


e , 

e* 

0 ^ 



e 

1" 

cr 

®11" 

^2 

CD 

ro 

ro 








P 

1 


radius of cylinder 

axial shortening per unit axial length 
see figure 7 


components of membrane strain 
= e‘/Ja 

= B^y'a = coefficients of Fourier series, equation (lO) 
defined by equation (,2k) 
see figure 5 

aspect ratio of buckles = i fl 

y/ X 

length of buckle in axial direction 

length of buckle in circumferential direction 

see figure 3 

number of circumferential waves 
see figure 3 

dimensionless critical loading = — - — 

E \fd 
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t 

U, V, V 


X 

y 

A, Cj D 
F, G, E, Q 
A', B', C’ 
Q* 

A", B",' C" 
Q" 

E 

E 

s 

E, F, G 
L, M, N 


K 

L 




S, S', R 


T' 

m m 

11 ^ 12 

T 

22 


thickness of cylinder wall 

con^jonents of displacement of a point on the middle sur- 
face of the shell 

coordinate in axial direction 

coordinate in circumferential direction 

coefficients, see equations (35) ^ (36) 

coefficients, see equations (37) ^ (38) 

coefficients, see equations (60), (6I) 

Young's modulus 

elastic constant of loading spring 

coefficients of the first fundamental form of the deflec- 
tion surface 

coefficients of the second fundamental form of the deflec- 
tion surface 

Gaussian curvature 

length of cylinder 

radii of curvatvire of deflection surface 
see figure 7 


membrane stress components 



k 
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a 


6 

6 * 

e* 

A 


total extensional and lending strain energy, respectively 


see equation (56) 



12(1 - y2)a2 


= a - h 
= p - a 

axial shortening of cylinder per unit length 


2jta 




= n 

= equation (k-3) 

Poisson’s ratio 

axial compressive stress 



dimensionless strain energy, see equation (52) 


dimensionless strain energy, see equation (57) 


X 


stress function for the membrane stress 
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EXPERIMEDTOAL EESULTS ON ELASTIC BUCKLING 


As for the 'buckling load of a cylindrical shell, the existing 
results appear sufficient for the present. But it will he useful from 
theoretical considerations to obtain experimental data concerning the 
difference between the so-called buckling load and the load after buck- 
ling, the characteristics of the buckling deformation, and the effect of 
the rigidity of the loading equipment etc. 

To clarify the post-buckling behavior, experiments were carried out 
with cellvilold cylinders. The restilts are shown in figure 1 for three 
kinds of rigidities of loading springs. It will be needless to say that 
the value of q. at buckling is much smaller than the value q. = 2 of 
the classical theory. For cylinders whose dimensions are shown in the 
figure, the post-buckling load is about one third of the buckling load. 

The clrciomferential number of lobes ranges from the value ^ to J, and 
their aspect ratio from O.65 to 1 . 0 . From the con^jarison of the three 
experiments, we can deduce that, with the Increase of the rigidity of the 
spring, (1) both the buckling and the post-buckling load Increase, 

(2) the value of k decreases, that is, the buckling wave form becomes 
more slender in the axial direction, and (3) the clrciimferential number 
of waves n increases. When the axial shortening is reversed, the load 
is kept nearly constant, and the value of k decreases, l.e., the wave 
length of the buckling lobe becomes longer, under the constant value of 
n. Such a change of the wave form will be proved theoretically reasonable 
later . 


CHARACTERISTICS OF BUCKLING DEFORMATION 


As for the buckling deformation of a cylindrical shell, there are 
two remarkable characteristics, other than those of the wave number and 
the aspect ratio of the lobes . The first is that the cylinder practi- 
cally always buckles locally, at least in comparatively long cylinders as 
used in our experiments. The longitudinal range over which such buckling 
occurs extends over I.5 times the axial wave length of the lobe; while 
circumferentially, it extends over the whole cylinder only under the con- 
dition of complete rotational symmetry of both the cylinder and the 
loading equipment. As will be considered later, the local buckling is 
attributed to the energy condition. 

The second characteristic is that the buckled surface, which is 
greatly displaced from the initial cylindrical surface, is nearly devel- 
opable. That is, the buckled surface is a regular repetition of the 
polygon ABCD (fig. 2 ), the ridges AB, BC, CD and DA and the trough 
BD being all nearly straight, and the surface ABD and BCD nearly 
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plane. Accordingly, the buckled svirface is approximately a concave 
polyhedron constructed by connecting nearly plane triangles such as ABD. 
Consequently, the Gauss curvature 


K = 


RlRg 


( 1 ) 


is vanishingly small throughout the deformed surface, where R^ and 

R^ denote the principal radii of curvature. Denoting the fundamental 

metrics of the first and the second kinds of the deformed middle surface 
by E, F, G and L, M, N, respectively, the compatibility condition 
is given by Gauss' characteristic equation 

1 2 5^g 

2 ^xSe ^q2 ^^2 

2 2 I 

H = EG - F J 

under the neglection of the second order terms regarding the derivatives 
of E, F, G. The left hand side of equation (2) is considered to be 
small, and as equation (2) is satisfied by the undeformed initial surface, 
on which the Gauss curvature vanishes completely, the expression of the 
left-hand side of equation (2) vanishes also for the initial s-urface. 
Consequently, when the deformed surface of the circular cylinder is nearly 
developable as in the case of buckling under axial compression, the defor- 
mation can be said to be almost inextensional. The deviation of the buck- 
led middle surface from the complete developable surface will depend on 
the thickness of the shell. For 'a limiting case where the thickness tends 
to zero, i.e., for an ideal membrane, the buckled surface is supposed to 
become e complete concave polyhedron. 

In case a completely developable surface is realized, the cross sec- 
tion through the point A in figure 2 will form a regular polygon as is 
shown in figure 5(a). Assuming the polygon to be n-slded, then the length 
of each side is 2rta/n, because the deformation is conpletely Inexten- 
sional. Hence, from the figure 


= LN - M = KH 


( 2 ) 


h = — cot — = a 
n n • 



and the displacement at the trough is given by 


6 = a - h = a 

,_2 


(5) 
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Since 

p = ~ cosec ^ = a/l + — — ^ 

^ " ■ V 6ny 


the displacement at the peak is 





Putting the aspect ratio k of a buckling lobe which has the wave 

lengths 21 and 21 in the circtimferentlal (y) and axial (x) 

y X 

directions by 


k = 


y 


( 5 ) 


we have 


^ _ ita 


(6) 


by z = jta/n. As we have 

y 


lx' = /zx^ - (6 + 6')^ 


From figure 3(b) ^ which shows the axial section of the deformed cylinder 
through A in figure 2, the axial shortening e of the cylinder per 
unit axial length is given by 


_ ^x - Zx' i 1 (6 + 5 * ) 




X 


* 2 12 
X 


and hence^ from equations (3), (^) and (6), 


e = 


-2 n 2 

jt k 

8 2 

n 


( 7 ) 


From equation ( 7 )^ i't is seen that the axial shortening of the cylindri- 
cal membrane, due to the geometrical deformation. Increases as the lobe 
becomes flatter in the axial direction and the circumferential wave num- 
ber decreases. Although the compressive load required for such a buck- 
ling into a conipletely developed surface vanishes, it does not vanish for 
an actioal shell owing to deviation of the actual deformation surface from 
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being con^Jletely developable because of its bending rigidity. The fact 
that a relation similar to equation (7) holds also for the shell, as will 
be shown later, is considered to be an evidence for the buckled surface 
to be approximately developable. 

The experimental fact above-mentioned, that the buckled surface of 
the cylindrical shell is approximately developable, is verified from the 
following theoretical consideration. If the higher order terms are 
neglected in the expression of L, M, N, with respect to derivatives of 
displacement components, it follows that 


L J ^xx ^ = ""xy ^ ^ + "'yy) (Q) 

where we write w, x, y instead of w/a, x/a, y/a, respectively. 

If, in the following calculation, one likewise puts u, v in place of 
u/a, v/a, to make all quantities dimensionless, then by equation (2) 



w 


XX 


(1 + Wyy) _ 



(9) 


To express the almost developable surface as shown in figure 2 by the 
displacement w, it is necessary to take a comparatively large number 
of terms of the Fourier series 

w = cos i?oc cos jpy (lO) 


where 


6 . . 

- = — 

ij a 


( 11 ) 


5^. denoting the deflection of each component of deformation, and 


7 . = 


2ita ^ _ 2jta _ ^ 

2Zx 2Zy 


( 12 ) 


n being the circumferential wave number. 

As a simple case of equation (lO), if we assume 
w = f^ + f]_ cos cos My ^2 


(15) 
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which Is nothing hut a wave pattern 

cos Toe cos (ly 

which has the nodal lines parallel to the x- and y-axes, superin 5 )osed 
by a wave pattern 

cos (Toe + |iy) cos (Ax - (jy) 

which has the nodal lines parallel to a line having an angle 
t arc tan ^ with the x-axis, then we have 

a K = A^ ^ 

2/1 P2 \ i222 

A ^ “ ^^2/ 2 ^1 ^ ^ 

. 2 2 2 

4f2 A M. COS 2Ax cos 2ny + 


%^2 


2 

A |i 


(^cos Ax cos 5^ly + cos 37vx cos ny) 


( 14 ) 


Only the first and second terms in equation (l4) can be made vanishing, 
and, if we choose the values of f^ and f^ so as to make these terms 

zero, it follows that 


= il 

1 n2 



(15) 


It is noteworthy that equation (15) has a form quite similar to 
equation (3). H Is assumed in the above derivation that the series (13) 
has only three terms, and a choice of f and f is made which makes 


K as small as possible. If we apply the infinite series (lO) generally, 
there exist the values 





= const 


(16) 
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•which make K equal to 0 completely. Therefore, the number of terms 
of the series must be increased, as the thickness of the shell decreases, 
and, for the limiting case of the ideal membrane, it will be necessary to 
take an infinite series. 

For a flat plate, the Gauss curvature K of the deflection surface 
cannot vanish unless a.1 1 the coefficients f^j vanish, by reason of non- 
existence of the term w in the right-hand side of equation (9) • This 

20c 

shows that no developable siirface exists in flat plates other than the 
initial plane, contrary to the case of the circular cylindrical shell. 
Consequently, when a flat plate performs finite displacement, severe 
bending must be accompanied with large extension of the middle surface. 
This is the reason why a flat plate and a cylindrical shell exhibit sub- 
stanially different behaviors diuring buckling. 


GENERAL BUCKLING 


Actual buckling is local in all cases, as was described in the last 
paragraph, except when a cylinder has a particular length. However, in 
order to treat the local buckling, it is first necessary to investigate 
the phenomenon of general buckling. 


Before proceeding to the subject, it is necessary to prove the exis- 
tence of the stress function even in the state of large deformation after 
buckling, because it has hitherto been proven only for small deformations 
The existence of the stress function has usually been considered to mean 
that the third and fo'urth terms can be neglected as compared with the 
first and second terms in the first and second equations of the equilib- 
rium equations (ref. 6) 


St St 

—11 + _ 

Sx Sy 

St St 

12 + 22 _ 

Sx Sy 

St St 

Sx Sy 


T L 


15 


- T M = 0 

23 a 


T 

13 a 




a 


T L + 2T 
11 


M + 

12 a 


T, 


22 





(17) 


Physically this concept means that bending can be neglected against exten- 
sion, and that the equations are satisfied for the usual extenslonal defor- 
mation. For the almost inextensional finite deformation, however, the 
existence of the stress function cannot be readily concluded in a similar 
manner, because T^^, T^^ '^22 small, and L, M, and N 
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large after buckling. Nevertheless, it can be derived by some other rea- 
soning when the deformed surface remains developable. 


Solving the first and second equations in (I7) with respect to 
and '^2^* have 



2 

For LN - M vanishes on the developable surface. Consequently we can 
conclude the approximate existence of the stress function for the buckled 
state of the circular cylindrical shells. 

In consequence, we can put 

( 20 ) 

using the stress function X. The compatibility equation can be written 

^ f A f s 

V X = - ^LN - M J (21) 

neglecting higher order infinitesimals . The deflection w must be chosen 
so as to minimize the right-hand side of equation (2I) . If we take equa- 
tion (15) which could nearly, thou^ not completely, satisfy this condi- 
tion we have, according to equation ( 21 ), 




" 7 ' 2 ~ "" 2\ 2 - ^ 1 ) ^ ^ 

+ n / 

? \ 1 A^ 2 A^q^fn fo ^ , 

— f \ cos 2 Ax - -= fn cos uy - ■ cos Ax cos 5 ^y - 


.2 2 . 


2 2 


/ - g— | vp 5 Ax cos py - i ^ - f/ cos 2Ax cos 2py - || y^ 

V9A + p j r / 


( 22 ) 
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where cr is positive in compression. Obtaining , ^12^ ^22^ f^om 

equation (20 ) , and then the strain of the middle surface e , e , e , 
^ " 11 12 22 

we have 


Ux = 


2 p 2 p 

- d 2i_f^ A f ^ 

~ E ~ a l~lv 2 


2 p ? 2 

=v| - ^ f - ^ f. + f 


8 1 




( 25 ) 


according to the relations 


1 2 

e = \iv + — w 
11 ^ 2 X 


e = V 

22 y 


2 

_ w + i W 

2 y 


and neglecting the periodic terms. 


If we put 


f 


1 




(2^) 


by means of equation (3) or ( 5 ), and g^ may be considered to be 

approximately constant for the post-buckling siarface. The concrete idea 
for the mechanism of buckling deformation will be obtainable later by 
Introducing these parameters. equation (5) 


k = A = 

n 


(25) 


and therefore the axial shortening of the cylinder per unit axial length 
is obtained from equation (23) as 




(26) 


where 


a = 


E fa 


e = 2 
E 


q = n^^^ 


a = 


12 




(27) 

(28) 
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And from the condition Vy = 0, ve obtain 




(29) 


Equation ( 26 ) shows that the total axial shortening of the cylinder 
consists of two parts, one of which is the elastic compressive strain 


e^ = q 


(50) 


caused hy the mean compressive stress a , when the shell is assumed not 
buckled, and the other is the shortening from the change of the geometri- 
cal shape of the siarface 


e 

2 


n 




(31) 


when the shell is assumed not to be subjected to the compressive stress. 
Similarly, equation ( 29 ) shows that the same circumstances axe also satis- 
fied for the mean radial displacement f^. The fact that equation (51) 

has the same form as equation ( 7 ) will be considered to be the consequence 
that the buckled surface is similar to the concave polyhedron mentioned 
previously (fig. 3(a)). The discrepancy in the value of the coefficient 


1 2 
8 ®1 



from that of 



will be a measure of the deviation of 


the surface from the conqiletely developable surface . 

Denoting the total elastic energy, the extensional, and bending 
energy by W, and W^, respectively, and putting 


W = 2rtaL — Etall 
2 

W = 2jtaL — Etall^ - 

e 2 ® 

W = 2jtaL i Etall-h 

b 2 “ 


(52) 


L being the length of the cylinder, E the Young's modulus, we ha.ve 


n = 


(33) 
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where, hy the use of equations (13) and ( 22 ), 


n = -^ + (^ — Y = i Q + 4^ 


\ e £/ tj 


n = ^ jjg ^2 


Q = Ag^^ + Bg^^g^^ + Cgg^ - Dg^^gg + Cg^^ + Fs^ 


+ 1 


B = k 


(k^ + 1)^ k(9k^ + 1)^ Yk^ + 9 )' 




k . 

" 2 \2 16 

.k + li 


" = l 


' 2 ^2 
vk + 1 j 


H = 2 (k^ + 1) 


The extensional energy consists of the elastic energy q due 

to the mean stress q and the energy — Q due to the deviation from the 

mean. The bending energy II^ does not depend on the circumferential 

wave number n. The cylindrical shell is distinguished from the flat 

2 

plate by the fact that the negative term -Dgj_ g2 exists in the expres- 
sion Q, and in consequence of this, Q, can diminish to any extent by 

increasing the number of terms, provided that the value of g. . is chosen 

J- J 

appropriately. Contrary to U^, the bending energy increases with 

the number of terms of the series representing w. Consequently, even for 
the finite deformation after buckling, the extensional and the bending 
energy have the same order of magnitude, and this also proves that the 
deformation is almost inextensional from the energy viewpoint. 
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Since the buckling load of circular cylindrical shells is to be 
determined from the energy level before and after buckling, the loading 
condition must be prescribed. We can consider constant end displacement, 
constant end load such as by dead weight, or other Intermediate states. 
Among these conditions, the fundamental case of constant end displacement 
will be treated first. In this case, the axial shortening e must be con- 
sidered as the independent variable. Therefore, from equation (26), the 
first equation of (3^) becomes 

H = Q« + - 2e ^ 

■ e 2 11 

n 


where 

Q = A‘g^^ + + C'gg^ - Dg]_^g2 + Cg^^ + Fg^^ (38) 


A 

A* = A + ^ 

64 

4 

B’ = B + 

64 

C* = C + i 4 - 

64 , 


Since the work of the external force is zero under the condition 
e = const., we have the equilibrium condition 


(59) 




6W = 0 


or, assuming the dimensions of the shell to be given. 


^ = 0 ^ = 0 (41) 

Sg^ 


under the condition that e, il(n) and k are constant, because U is 
considered to be 


(iv2) 


n = n(e, n(n), k, gjL, gg) 
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From equation (4l) we obtain 

i + 2B»g^^p^ - 2Dg^p + 2cJ + 2G - I e ^ = 0 

i r2B*g,^p + l<.C*gT^p^ - Dg^ + 2Fpl + 2Hp - e p = 0 

t. 2 L 1 1 1 J n 


(^ 5 ) 


82/^1 - P 

Eliminating e from equation (43), we have an equation of the third 
order with respect to p, namely. 


4 (b* - C‘)g^^p^ - 4Dg^p^ + 

r(aA‘ - 2B‘)g^^ + 4C - 2F + (4G - 2H)] p + Dg^ = 0 (44) 

Substituting the value of p, obtained from equation (44) with parameters 
k, T) and g^^, into the first equation of (43), we obtain e, and then 

q from equation (26) . The relationship between q and g^ is shown in 

figure 4 for each case of k=0.2, 0.6, and 1.0, Though the value of g^ 

for the post-buckling equilibrium state is not determined from these 
curves^ the fact that all curves lie within a definite range of g^ and 

show a similar shape, indicates the reasonability of the assumption (24). 

The relation between q and e is shown in figure 5 with k and 
T] as parameters. The minim'um value of q, corresponding to the branch 
points of the curves exhibiting the post-buckling state, is the critical 
load by the classical theory 


q = 2 


as is well known. 

It is important that the curves in figure 5 reveal the condition 
H = min. (only mder the restriction k) , ii = const, but do not give 
the minimum condition of II ^e, ri, k, g^, g^^ with respect to each 

of its independent variables. Therefore, the inference that the smallest 
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value of q in figure 5 will "be realized after "buckling, is erroneovis. 

In other words, the fact that q is negative and has a large absolute 
value for small values of k means that the geometrical shortening ( 51 ) 
is greater than the total shortening e, and therefore the extenslonal 
- energy 11^ is Increased tinder such conditions, and II does not become 
minimum. 


As was stated above, for the determination of the post-buckling 
state, it is necessary to determine the values of e, tj, k, g^, and 

gg, and hence of q, corresponding to the minimum value of the hyper sur- 
face ^e, T], k, g^, g^^ with regard to the variables e, T, k, 

g^ and gg. For this purpose, we need to inquire into the minimum value, 

with respect to t] and k, of the minimum of II in a section tj, 
k = const. Figure 6 shows the value of the total strain energy II cal- 
culated from figures ^4- and 5^ taking k and as parameters. Hence, 
this calculation gives the minimum of II for each section. Each curve 
in figure 6 has the form as shown in figure TC^). For e < e* in the 
figure, buckling is Impossible, but for e*< e< e^, it is not impos- 
sible if some energy is supplied from the exterior. It is for e > e^ 

that buckling takes place in general, but for this a certain amount of 
energy barrier must be jumped over. At e = ej_ > e^, for example, the 

energy barrier corresponding to CC must be jutrped over in order that 
the equilibrium position shifts from C to C" by buckling. This pro- 
cess can be illustrated as in figures 7(l3) and (c) with g-j^ or q as 

the Independent variable. As the buckling load corresponding to e^ 

Increases, the energy barrier to be jumped over decreases, and the energy 
difference before and after buckling increases. 


From figure 6, it is shown that the value of g^ and g^ corre- 
sponding to 6q Is almost constant. Independent of k and q; this is 

worth noting as indicating that the starting assun^ition (24) is reason- 
able. 


Denoting the minima of e^ with respect to 

k in figure 6 by e , and the values of II 

o,m-' 

of q at S (refer to fig. 7 ) corresponding to e 


T) for each value of 
at S and S' and that 
■by 


n n- 

o,m “ o,m^ o,m^ 
and T]^, respectively, we can obtain the manner in which these quantities 

change against k (fig. 8) . The value of ^ seems to show a mini- 
mum II near k = 0.2, and the corresponding value of e is 

o,m,m ' JT o 

about 1.3 and that of 0.02. Therefore the buckling load in this 
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case is q = 1.3, while the load after buckling is estimated roughly as 
q = 0.9« The energy barrier to be jumped over is approximately 

n = 0.l4, which is equal to about one tenth of the energy just before 
b 

buckling. Generally, buckling occurs for q > 1.3, since the value 

q = 1.3 corresponding to m m miTilminTi buckling load. The 

greater the buckling load, the smaller the energy barrier, i.e., the more 
likely buckling will occur, while the greater the energy difference before 
and after buckling, i.e., the greater the sound or vibration is generated 
with buckling. The problem, at what load over the minimum buckling load 
the actual buckling will occur, is to be decided by the Initial deforma- 
tion of the shell and the conditions of the experiments. With a constant 
value of e greater than 1.3, the value of k and consequently of tj, 
corresponding to the minimum of the energy is somewhat greater than that 
for the case e = 1.3. Therefore, as the actual buckling load increases 
over the minimum buckling load, the lobes have a tendency to become short 
and thick in the axial direction, and the circumferential wave number 
Increases. The energy curve for any constant e is very flat in the 
neighborhood of its minimum, hence the post-buckling state cannot be deter- 
mined sharply, but will be scattered over a certain range. At any rate, 
the lobes in the case of the general buckling aire found to be very slender 
in the axial direction. 

Though the general buckling was not realized in experiments with 
comparatively long cylinders, it will be realized if a cylinder of a par- 
ticular length, corresponding to 1.5 times of the wave length, is used. 


LOCAL BUCKLING 


According to the experimental results, actual buckling is local, 
consequently there exist two parts, one of which is buckled and the other 
unbuckled. The occurrence of local buckling is explained qualitatively 
from the following energy consideration. If local buckling were assumed 
to take place at the minimum buckling load of general buckling, the load 
supported by the buckled part would decrease abruptly; and in consequence 
of the fall of the load, the elastic energy of the unbuckled part will 
decrease more than the energy Increase of the buckled part. Hence the 
total elastic energy would become smaller after buckling under fixed end 
displacement. This means that the minimum buckling load in the case of 
local buckling will be found below the general buckling load. The occur- 
rence of local buckling, on the other hand, is Inferred from the genetic 
viewpoint that the stress-strain relation after buckling (fig. 5) is such 
that the velocity of propagation of the buckling deformation is nearly 
zero. 
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In this paragraph we shall consider this problem from the energy 
standpoint on the basis of the results on general buckling. Dividing the 
entire cylinder into the buckled and the imbuckled parts and writing 


Buckled 

part 


Unbuckled 

part 


Length 

Mean compressive stress 

Compression 

Elastic energy 




we have 


Entire 

cylinder 

L 

a 

e» 

W 


L ■— ^1 ^2 

= q = Ee^- 

W* = W^* + Wg’ 


Le* 


Vi 


' + L e ' 
2 2 


(^ 5 ) 

(46) 

(47) 

(48) 


where 


q. = 

Eliminating e^' and e^’ 




from equations 


( 14 - 6 ), (26), and ( 48 ), 


(49) 


we obtain 


= q + 


_2 ^ /I 
L n V8 


. 2 ^ 1 2> 
=1 4 ®2 > 


(50) 


Equation (50) shows that the total axial shortening depends on the 
length L2 of the buckled part. As we can put 

Lg = 1.5(2Zx) 


(51) 


20 


MCA TM 1390 


from the experimental results, we have 


T 3ita 
^2 - ^ 


(52) 


from equations (5) and (l2) . Therefore equation (50) is written as 


e’ = q + 


3rta k 
Lnk T 


i i g 


(53) 


The second term of the right-hand side of equation (53) represents the 
geometrical shortening due to local buckling. 


The total elastic energy is the sum of the elastic energy of the 
unbuckled part . 


= 2jratLi I Ee^' = 2jrali ^ aq^ 


i^) 


and that of the buckled part 


Wo * = W 


(55) 


where W is given by equations (52) and (3^) . Hence, putting 


W' = 23taL I Etall 
W^' = 2itaL I BtaU' 
W 2 * = 2jtaL i Etanp* 


2 J 


(56) 


it follows 




Lp r 

n ' = 

2 L 


0, + q^ + 


" ^2) 


or by means of equation (52) 


n- . n^' -H = g| 


i 9" H- (og," + Hg/) 


+ q. 


(57) 


(58) 
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Considering e* to be the independent variable, eqxiation ( 58 ) is written 


as 


H' = 


3ita 

Lnk 


i «” + + Hg^^) 

" (I ^ t 


gi2 _ 2e* 

Lnk 11 


by equation ( 53 ), where 


Q" = A"g^^ + B"gi^g2^ + C"g2^ - Dg^^gg + Cg3_^ + Fg2^ 


( 59 ) 


(60) 


A" = A + 


3jfa 

Lnk 6k 


B" = B + ^ 

Lnk 6k 

C" = C + ^ 

Lnk 64 


( 61 ) 


The minimum condition for W, or = 0 and = 0 under the 




dition e', k = const, is given by 




con- 




^ + 2B"g^^p^ - 2Dg^p + 2C 


+ 2G - ^ = 0 

2 Ti 


i r 2 B"g 2p + 4C"g^^p^ - Dg^ + 2Fp 


Tj^ L 


"1 • - -1 - ---i 

Eliminating e' from equation (62), we have 

23 


+ 2Hp - e' p = 0 


(62) 


4(B" - C")g 3 _ - tog^P + 

[(aA- - 2B")g^^ + 4C - 2 F + t]2 (4g _ 2H)J p + Dg^ = 0 ( 63 ) 
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which is equivalent to equation (44), 


as the relation 


B" - C" = B* - C‘ = B - C 
8k" - 2B" = 8A* - 2B' = 8A - 2B 


(a) 


is satisfied hy means of equations (61) and (39) • Consequently, we can 
use the same results regarding g^, p, and t) as were calculated in the 

case of general buckling. Comparing e’, obtained from the first equa- 
tion of (62) by means of calculated results, with e derived from the 
first equation of (43), we have 


e* 



(65) 


This equation reduces to equation (55) by virtue of equation (26), so 
that the value of q in the case of local buckling is the same as that 
of q in the case of general buckling. Hence, we need not carry out the 
calculation for the local buckling under the condition e* = const, all 
over again, but only calculate e* and H‘, which satisfy the minimum 
condition of energy with respect to ^ and by equations (53) and 

(58), respectively, using the results calculated for the case of general 
buckling. Local buckling differs from the case of general buckling in 
that it depends on the dimensions a, L, t and the Poisson's ratio v 
of the circular cylinder, and therefore on the value of n. 


Using the values a = 111 mm, L = 460 mm, t = 0.4 mm, and v = 0.4 
for the shells used in the experiments, the e' ~ q and the e' ~ II’ 
relations are calculated as shown in figure 9 and figure 10, respectively. 
As in the case of general buckling, we can define for local buckling e^'. 


H H-’ 

o ^ o 

mum of II ’ 
o 


and e ' , 

o,m^ 

with fixed 


H’ , H'- n corresponding to the mini- 

0, m"^ o,m'^ m ° 

k. The manner in which these quantities change 


against k is shown in figure 11. The minimum e’ or H’ of 

^ o,m,m o,m,m 

e' or H’ exists near k = 0.7, and e’ = 1.04, and the 

o,m o,m ’ o,m,m ’ 

corresponding value of q is approximately 0.10 (or n = 9»40). The 

post-buckling load for e' = ®'q = 1*04, k = 0.7 and tj = 0.10 is 

obtained as about q = 0.7, and the energy barrier H’- - II' 

^ ^ o,m,m o,m,m 

about 0.07. The post-buckling state when buckling occurs at any load 

above the minimum buckling load q = e' = 1.04, will be determined 

^ q,m,m ’ 
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by the minimum energy condition at e* corresponding to the buckling 
load q. This state, however, does not differ greatly from that of mini- 
mum energy represented in figure 11. 


Comparing figure 11 with figure 8, it is seen that the minimum 


energy II' 


o,m,m 


for the local buckling is much smaller than the mlnl- 
for the general buckling, and this means that the 

curve shows a 


mum energy II 

^ o,m,^ 

local buckling is more liable to occur. As the H' 


o.m 


sharper minimum than that for H 


o,m^ 


the post buckling state for local 


buckling will not be scattered as widely as in general buckling. 


The above calculation was carried out under the asstanption that the 
end shortening is constant before and after buckling. Hence the results 
may be compared with expierlment 1 (fig. l) perfonned with a spring of 
large rigidity (table l) . The agreement of the value of k is good, but 
that of the buckling load, the load after buckling and the circumferential 
wave number, is not very good. 


LOCAL BUCKLING WITH ELASTIC LOAD EQUIPMEnn' 


Actually, it is impossible to realize the condition e or 
e' = const, in buckling by usual testing machines. The absorption of 
energy by the testing machine is inevitable due to the finite rigidity of 
the structiural material of the testing machine, even when springs or oil 
pressure is not used. The absorbed energy increases with the decrease of 
the rigidity of the machine, and it is natural that the interchange of 
energy between the machine and test specimens have some effect on the 
buckling phenomenon. In the present paragraph, a case in which the rigid- 
ities of both the loading equipment and the testing machine are replaced 
by a spring, will be analyzed and compared with the experimental result. 

It is assxamed that a shell undergoes local buckling as in the pre- 
ceding paragraph. Denoting the shortening of the cylinder per imit length 
by e'(e' = e*/\^) as before, the length, the unit con^jression and the 
rigidity of the spring by L , e (e = stn<i Eg, respectively, 

and the length and the unit compression of the total system containing 
both the cylinder and the spring by L-fc and 

respectively, we have 


L^ L + Lg 


( 66 ) 


L,e, = Le* + 

L e 

t t 

s s 

e = at 


s 

E. 


(67) 


( 68 ) 
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Substituting equations ( 53 ) and ( 68 ) into equation (67), the shortening of 
the total system is given hy 


e. = (1 + ^ + 

L t \ L Eg / 


3Jta ^ fl 2 

Lnh n V 8 1 



As the elastic energy of the spring is given hy 


(69) 


w 


S 


= I ^S^S^ 


2 

s 


or 

W = 2 itaL i Etallc, 
s 2 ^ 

H ^ 2rtatE _2 
s “ L Eg ^ 


(70) 


the elastic energy of the total system is given hy 


¥" 

n" 


2 jtaL i Etall" 
2 


37 ta 

Lhk: 


Q + 




+ 


^s 2 rtatE \_2 
s / 


(71) 


using equations (36) and (58). The energy carried hy the spring is pro- 
portional to L and 


The buckling in this case is to he considered with e^ as an inde- 
pendent variable, which was e' in the case of local buckling. Substi- 
tuting, therefore, q of equation (69) into equation (71), we obtain 


H" 


3 Jta 

Lnk 


^ Q' " + (dgp^ + Hgg^) 


+ 


1 

1 + (LsAX2^atE/Eg 



2 



't Lnk 11 



+ 


1 

if 



(72) 
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Q 


I t t 


= A* • + B' • + C« « - Dg^^g^ + Cg^^ + Fgg^ (T3) 


B 


t » * — 


^ 3jra/(Lnlc) ^ 

1 + ^L^/L)^2«atE/Egj 64 

g ^ 3jta/ (Lnk) 

1 + (l^/l )''2jtatE/Eg; 64 


= C + 


3rta/(Lnk) 


1 + (L^/L)(2rtatE/E^) 64 


(T4) 


The equilibrium condition under constant n and k is given by 


5II" 


= 0 


5it" 

Sgo 


= 0 


or 


4A’ + 2B' ' 'g^^p - 2Dg^p + 2C 


_1^ 

_2 


2B 


2 2jtatE 4 

1 + — 2- 

L E= 


« ' I rr 


gj^P + 4C ’ 'g^^p^ - Dg^ + 2Fp 


(V9 


L 2rtatE 4 

1 + -^ 

L E, 


e 2 

p + 2Hp = 0, 


;■ (75) 


from which, eliminating e^, the same equation of the third order with 

respect to p as equations (44) and ( 63 ) is derived by means of equa- 
tion (?4) . Hence, the calculated results of q., 4 , k, and p in the 
case of general buckling may be used also in this case. It is convenient 
to use equation ( 71 ) instead of equation ( 72 ) for the calculation of 
energy, and equation ( 69 ) in place of equation (75) for the calculation 
of e^. Figures 12 and 13 show the relations H" and 
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n" ~ e’, respectively, for a spring with. Lg = 46 mm, = 5^680 kg 

used in experiment 2 and for the same cylindrical shell as in the last 

paragraph, while figures l4 and I 5 show those for a spring with 

L =43 mm, E = l,l60 kg (experiment 5*) The minimum buckling load 
s s 

is given by the value of e' or q corresponding to the minimum 

with respect to k and ti of e, , which is the value of e, 
t,o,m,m t,o^ t 

corresponding to the intersecting point of the post-buckling energy c\irve 

with the parabolic curve in figures 12 and l4. 

As e is constant before and after buckling in o\ir experiments, 

the buckling load and the post-buckling state will be determined by con- 
sidering these energy charts \jnder the condition e^ = const. For any 

value of e, above e, in figures 12 or l4, the corresponding 

t t,o,m,m JT a 

value of H" is obtained on the parabola, hence the value of e’ or 

the buckling load q is obtained from figures 13 or I 5 . The post- 

buckling state is determined by the following process. The values of k, 

T](n) and H" which give the minimum energy for the given e^ are 

determined from figures 12 or l4, then the corresponding value of e’ is 
determined from figures 13 or I 5 , and finally, the post-buckling load q 
corresponding to these values of k, t] and e’ is determined from 
figure 9. 

The results calc-ulated by the above procedure for the shell and the 
spring in experiments 2 and 3 are shown in table 1. In comparing with 
experimental results, the qualitative agreement between the theory and 
the experiment seems to be satisfactory. The reason that the observed 
minimum buckling loads in experiments 2 and 3 are not so small conpared 
with those in experiment 1 may be due to the magnitude of the energy 
barrier. Practically, the supposition that the buckling in experiments 
2 and 3 might occur at greater difference of the energy before dnd after 
buckling compared with experiment 1, or at energy barrier closer to 
experiment 1, is confirmed by the fact that the sound on buckling was 
greater in the former cases than in the latter case. It may be presumed, 
consequently, that the occurrence of buckling is governed both by the 
magnitude of the minimum buckling load and the energy barrier. 


deformation and stress distribution 


In order to clarify at what part of the buckled shell the conpres - 
slve load is supported, we calculated the stress distribution (fig. I 7 ) 
for the deflection w due to equation (13) (fig. I 6 ) . According to the 
result obtained for the post-buckling state q = O.7, k = O.7, and 
T = 0.10 (n = 9«4), the compressive load applied to the cylinder is found 
to be supported by the cross-hatched part. 
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CONCLUSIONS 


The huckled sijrface of a circiolar cylindrical shell under axial 
compression is considered to approximate a particular developable surface, 
so that the deformation in this case is almost inextensional and finite. 
The mechanism of buckling of a circular cylindrical shell can be said, in 
effect, to be a phenomenon of a kind of "Durchschlag." 

When the cylinder is loaded by a specified end shortening, the state 
realized after buckling is determined by the minimum condition of the 
hypersurface representing the elastic energy with respect to a1 1 variables. 
In consequence, there exists the minimum buckling load which is determined 
by the equality of the energy before and after buckling and is smaller 
than the critical load by the classical theory. Actiial buckling takes 
place above this minimum buckling load, by jun^iing over a certain magni- 
tude of energy barrier, which decreases with the increase of the buckling 
load. The buckling load will be considered to be governed by both the 
rainimimi buckling load and the energy barrier. 

From the above viewpoint, local buckling is more likely to occur than 
general buckling, the minimum buckling load and the energy barrier for the 
former being smaller than those for the latter. In the case of local 
buckling, the minimum buckling load decreases as the rigidity of the 
loading spring decreases. Similarly, the fact that the minimum buckling 
load for the local buckling is smaller than that for the general buckling 
is attributable to the spring action of the unbuckled part. 

Finally, it will be remarked that scale effects exist for the local 
buckling as is seen in equations (75) and ( 75 ), l.e., that the buckling 
strength is Influenced by the absolute value of the length, radius, and 
thickness of the cylindrical shell even under the condition of their con- 
stant ratios. A paper concerning this effect will be published soon 
(ref. 7). 

The author wishes to express his appreciation to Professor M. Yoshiki 
for discussions. This research was carried out by the financial aid of 
the Ministry of Education. 
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TABLE I 


. 

General 

buckling 

. 

Local buckling 

Loading with 
constant end 
displacement 

Loading with spring 

Exp. 1 

Exp. 2 

Exp. 5 

Minimum "buckling 

Calculated 

1.3 

1 . 0 k 

0.80 

0.80 

load (q) 


Observed 


0.8 

0.76 

0.72 

Energy "barrier (U, 

n', n") 

O.lL 

0.07 

0.20 

0.10 ' 


a 

Calculated 

0.9 

0.7 

0.55 

0.70 



Observed 


0.4 

0.24 

0.22 



Calculated 

0.2 

0.7 

0.75 

0,80 

Post-buckling 


Observed 


0.64 ~ 0.68 

0.93 

0.95 

state 

T](n) 

Calculated 

0.02 

0.10 ( 9 . 4 ) 

0.085 (8.6) 

0.085 (8.6) 



Observed 


6 ~ 7 

6 

5 



Calculated 



1.0 

1.25 



Observed 



3.2 

4.4 


ru 
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Figxire 1.- Experimental results 
aspect ratio k of the bucklini 

being g--^, e' 
t = the thickness of the shells, 
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Figure 4.- Relation between the stress q and the deflection g^ after 

buckling. 
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Figure 5.- Relation between the stress q and the unit end shortening 

e for general buckling. 
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0.2 


0.1 


0 


Figure 8.- Relation between the minimiirn II q of the elastic energy, 

which is equal before and after buckling, with respect to t) and the 
aspect ratio k of the buckling lobes, for general buckling. 
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Figure 9.- Relation between the compressive stress q anij unit end 
shortening e' for local buckling. 














0 


i 


2 



3 


Figure 10.- Relation between the elastic energy n' corresponding 
to its minimum condition with constant k and and unit end 
shortening e* for local buckling. 
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of the elastic 


Figure 11.- Relation between the miriimuin H ^ 

energy, which is equal before and after buckUng, with respect to 
n and the aspect ratio k of the buckling lobes, for local buckling. 



2 


Figure 12.- Relation between the total elastic energy II” corresponding 
to its minimum condition with constant k and tj and the unit end 
shortening e^ for local buckling of the system consisting of the 
shell and loading spring, the constants of which are equal to the 
case of experiment 2. 









ro 



Figure 13.- Relation between the elastic energy II'' and the unit 
compressive strain e' only of the shell for the case of figure 12. 
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Figure 14.- Relation between the total elastic energy II" corresponding 
to its minimum condition with constant k and n and the unit end 
shortening e^- for local buckling of the system consisting of the 
shell and loading spring, the constants of which are equal to the 
case of experiment 3. 
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Figure 17.- Distribution of the quantity ^ being the 


axial stress, negative value of which represents compression. 
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